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Abstract. The fidelity metric has recently been proposed as a useful and elegant 
approach to identify and characterize both quantum and classical phase transitions. 
We study this metric on the manifold of thermal states for the Lipkin-Meshkov- 
Glick (LMG) model. For the isotropic LMG model, we find that the metric reduces 
to a Fisher-Rao metric, reflecting an underlying classical probability distribution. 
Furthermore, this metric can be expressed in terms of derivatives of the free energy, 
indicating a relation to Ruppeiner geometry. This allows us to obtain exact expressions 
for the (suitably rescaled) metric in the thermodynamic limit. The phase transition 
of the isotropic LMG model is signalled by a degeneracy of this (improper) metric in 
the paramagnetic phase. Due to the integrability of the isotropic LMG model, ground 
state level crossings occur, leading to an ill-defined fidelity metric at zero temperature. 



PACS numbers: 05.70.Fh, 02.40.Ky, 64.70.Tg, 75.10.Jm 
1. Introduction 

Beginning with ground-state overlap studies of Zanardi and Paunkovic pQ, the fidelity 
of quantum states has recently been used for investigating classical as well as quantum 
critical behaviour in various systems. The motivating idea behind this approach is 
simple, yet extremely plausible: The properties of different macroscopic phases of matter 
should be encoded in the structure of rather distinct quantum states. Hence, a suitable 
metric that can quantify how "different" two given quantum states are should be able 
to capture some signature of a phase transition (see |2j for a recent review of these and 
related ideas). 

The appeal of this approach lies in the fact that it is related to geometric structures 
inherent to the state space of the given quantum system itself. This was already pointed 
out in [3] for the case of pure quantum states, and a generalization of fidelity to finite 
temperatures was discussed in Fidelity itself and the corresponding geometric 
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quantities might thus serve as "universal order parameters" that reveal signatures of 
criticality at zero as well as finite temperatures. A related approach, proposing the 
use of the so-called fidelity susceptibility in order to identify and characterize quantum 
phase transitions, has recently been put forward by You et al. [5]. 

A further interesting feature of the fidelity-metric approach lies in the fact that 
it also applies to non-standard (quantum) phase transitions, like topologically ordered 
phases [H]. For such transitions, no symmetry breaking principles are at work, and no 
local order parameter can be defined. For current results on fidelity and fidelity metric 
approaches to topological order, see [HIH]. 

In the present article, we study the phase transition of the Lipkin-Meshkov-Glick 
(LMG) model within the fidelity-metric approach [9]. Originally, this model was 
proposed to describe excitations in simple atomic nuclei. In its spin-1/2 representation, 
it can be regarded as a quantum XY model with infinite-ranged ferromagnetic exchange 
interactions, where every spin is subject to an external transverse magnetic field h. This 
model shows a continuous phase transition from a symmetry-breaking, ferromagnetically 
ordered phase to a phase that is spin-polarized for zero temperature and high fields and 
crosses over continuously to a paramagnet at zero field and high temperature. We 
mostly study the isotropic case, being rotationally symmetric in the (x, ?/)-plane. This 
case is somewhat special due to the fact that the Hamiltonian consists of mutually 
commuting terms, and no "competition" between noncommuting terms (regarding e.g. 
symmetry) can arise. Our aim is then to obtain the Riemannian metric tensor field 
related to fidelity, defined on the model's quantum state space. As expected, we find in 
this metric a signature of the phase transition. The peculiarities of the isotropic LMG 
model lead to a number of remarkable properties of the metric: First, as a consequence 
of exact ground-state level crossings, the metric is not well defined on the ground state 
manifold, i.e., at zero temperature. Second, for finite temperatures, we find a very 
pronounced signature at the phase boundary, with a well-defined Riemannian metric for 
the ferromagnetic phase, and a degenerate tensor field (not being a proper Riemannian 
metric) for the paramagnetic phase. Third, the metric components can be expressed 
entirely in terms of derivatives of the free energy, suggesting a close relation to Ruppeiner 
geometry [lOj. These features should disappear for the anisotropic case, i.e., as soon as 
a noncommuting term is added to the Hamiltonian. 

Studies of the phase transition of the LMG model within the framework of fidelity, 
fidelity susceptibility and related concepts have been reported previously. In [TT], fidelity 
was used basically as an alternative means to obtain the phase diagram, whereas in [12] 
the fidelity susceptibility and its scaling behaviour were studied (see also [13] and [H] for 
related work). Yet, to our knowledge, the explicit calculation of the associated metric 
tensor field at finite temperatures is novel. 

The article is structured as follows: In sections |2] and |3] we give an overview of 
quantum state space and its underlying geometric structures. This will lead us to the 
concepts of Fubini-Study geometry in the case of pure states and Bures geometry for 
mixed states. In section |4] we introduce the isotropic LMG model, its simple solution 
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in terms of angular momentum states, and its exact thermodynamic solution, section 
[5] is devoted to the computation of the fidelity metric induced on thermal submanifolds 
and the Ricci scalar. The Fubini-Study limit is discussed in section |6} and remarks on 
the anisotropic LMG model can be found in section [7| A discussion of the results and 
an outlook on future work is given in section |8j 

2. Fubini-Study Geometry on Quantum State Space V{Jif) 

In this section we introduce a Riemannian metric on quantum state space which serves 
as a measure of distinguishability of quantum states. As a first step, following [T5] . 
we will introduce quantum state space as a base manifold of a certain fiber bundle. 
Then there exists a very natural (from a mathematical point of view) way to derive a 
metric on quantum state space from the scalar product on Hilbert space. Remarkably, 
this metric has an information-geometric interpretation, rendering it a useful measure 
of distinguishability of quantum states. 

Consider a quantum system defined on a Hilbert space J^. We denote by 

= {|^) G =^ I = i}cj^ (1) 

the subset of normalized Hilbert space vectors. Then it is well-known that the 
relevant physical information is contained in the transition probabilities | (t/'Iv?) p, where 
{if) G S{Jif). However, S{J^) contains redundant state vectors, and therefore is not 
what we would like to call the quantum state space: For a phase-shifted Hilbert-space 
vector 

m^e''\^), (2) 

it is obvious that K'?/''!;/))!^ = K-i/^ly?)!^, and lip) and \ip') cannot be distinguished by 
measuring expectation values of any observable acting on alone. Putting it less 
mundane, the invariance of transition probabilities under these U{1) transformations 
induces an equivalence relation \ip) ~ \ip') on S{J^). We denote by [ip] G V{Jif) 
the corresponding equivalence classes, where the projective Hilbert space V{J^) is the 
space of equivalence classes. The projective Hilbert space now defines our first version 
of a quantum state space. Note that, for finite-dimensional Hilbert spaces = C^, 
the projective Hilbert space V{J^) = CP^~^ is a complex projective space, which is 
well-studied in geometry. The projection mapping 

n:S{J^)^V{J^), (3) 

allows for a fiber bundle interpretation: S{,y^) V{J^) is a principal fiber bundle 
with structure group U{1) and quantum state space V{Jif) as its base space. The fibers 
7r~^ (['?/']) are one-dimensional subspaces of and are themselves isomorphic to U{1). 
Note, that V{Jif) is isomorphic to the space of one-dimensional projectors of the form 
|^)(^|. 

The Hilbert space possesses a geometric structure that originates from its scalar 
product in a straightforward way, 

{^\^)^G{i;,^) + in{^,^). (4) 
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Here, G{iIj, (p) and VL^ip, ip) are defined as the real, respectively imaginary, part of 
G : ^ X ^ R is a bilinear, non-degenerate and symmetric map. Due to linearity 
of the Hilbert space ^ , we can identify its tangent space T^M' at a point with ^ 
itself, T^M' = Jif. Thus, G can also be seen as a mapping from TJif x TJ^ to the 
reals, and indeed defines a Riemannian structure on the Hilbert space J^. Similarly, 
Q : Jif X —y R defines a symplectic form on Jif, and together with G it endows 
Jif with a Kahlerian structure. The interested reader can find more information on 
these geometric structures and their implications for quantum mechanics in references 
[T5| [16] . In the present article, we will be concerned exclusively with the properties of 
the Riemannian metric G. 

Our next aim is to carry over the Riemannian structure from S{J^) to the quantum 
state space V{J^). Clearly, the Riemannian metric defined in ^ is not invariant 
under the U{1) phase rotation But the metric structure in the projective space of 
equivalence classes cannot depend on these phases and should be defined accordingly. 
Here the bundle structure S{J^) V{J^) comes in handy. A connection on a fiber 
bundle introduces the notions of vertical (g V^) and horizontal vectors (g H^). Vertical 
vectors "point along" the fiber direction and are elements of the tangent spaces to the 
points in tt^-*^ ([■?/;]). So given a curve (— e, +e) 9 t i— > e'^*^*^ \ip), 9{0) = along the 
fiber 7r~^([-?/']), the tangent vector ^(e'^'-*-' \i'))\t=o = i^(0) \ip) at the point \'ip) spans the 
vertical vector space = iR. To describe the connection in terms of a 1-form, we can 
naturally make use of the Hilbert-space scalar product. Take as this natural connection 
{ipl') '■ T^S{J^f) — > C. Then the horizontal tangent space at a point \ip) G S{J^) is 
given by those vectors which are mapped to zero by this connection, 

H^^{\^)eT^Si.^)M^) = o}. (5) 

This is precisely the orthogonal complement to {ip) in S{Jif), yielding the decomposition 
T^^(^) = H.^®V^ (6) 
of tangent spaces of S{je). An element \ip^) G can now be written as 

\^'') = \^)-{^p\^m. (7) 

This enables us to define a bilinear mapping 

(■|-)[^]:%P(^)xT[^,]P(^)^C (8) 
on the tangent spaces T[^]'P(^) of V{J^) as 

(Pi|P2)m = (^ly'f)- (9) 

Here, \(Pi), Iv^^) are vectors which are pushed forward to Pi,P2 G T[^]P(^) by the 
tangent projection 71^, : T^S{J^) — > T[^]V{Jif) that gives the tangent vectors to the 
projected curves 7i(\Lp(t))) in the base space V{J^). Eq. ^ can also be written as 

{Pi\P2)m = {vi\y^2)-{vi\i^){i^\y^2). (10) 

This object is often referred to as the quantum geometric tensor [17j. By construction, 
it is invariant under the U{1) transformation introduced above. 
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Taking the real part on both sides of equation (10), we can now define 

g{P„P,)^^m\P,)[^]} = G{v?,^^). (11) 

as a Riemannian metric on the projective Hilbert space. For exphcit calculations, it 
proves useful to rewrite g by employing a local section V{Jif) — * S{Jif), [ip] i— > \ip)- 
This section induces a push forward 

%P(^) r^^(^), P ^ |d^(P)), (12) 

which allows us to write \(p) = \dilj{P)). Finally we obtain 

g{PuP2) = ^{{dij{P^)\dij{P,)) - (d^(Pi)|^)(^|d^(P2))}. (13) 

This Riemannian metric is usually called Fubini-Study metric. Remarkably, one finds 
that the distance corresponding to this metric is a distance in the information-geometric 
sense, telling how 'difficult' it is to distinguish between certain states by means of ideal 
measurements (see |18j for details). This metric is useful for studying quantum phase 
transitions at zero temperature, where only pure quantum states need to be considered. 



3. Bures Geometry on Quantum State Space Ai 

For the study of thermal phase transitions, we have to extend the fidelity metric to 
mixed states, i.e., to the space of density operators. This formalism is mostly due to 
Uhlmann [12]; the presentation in this paper mainly follows reference [T3] . 

Let Ai denote the set of density operators, which defines our second version of 
a quantum state space. First note that V{J^) C Ai, since we have an isomorphism 
[ip] I— s> \4'){ip\ for pure states. For mixed states, one can identify a fiber bundle structure 
by observing that any density operator can be written p = WW'^ . The so-called 
purification W G S{,j^^^) is an element of 

S{.^^^) = {We,^^^\\\W\\ns = l}, (14) 

the Hilbert-Schmidt space of bounded operators W : —>■ with unit norm 
||iy||Hs = a/ {W, W)iis that is derived from the scalar product 

{Wi,W2)ns = ti WlW^. (15) 

Now, what is this construction good for, and where is the bundle? The purification of a 
given density operator is not unique since, if W : defines a purification, 

then WV with V G U{.y^) purifies p as well. Here, U{,y^) denotes the group of 
unitary operators acting on the Hilbert space J^. So instead of considering equivalence 
classes over just the pure-state phases U{1), we now introduce a projection mapping 
TT : S{J^^^^) ^MhyW ^ p = WWl 

There is still a slightly technical obstruction to obtaining a well defined U (J^) 
bundle. For general density operators, the "fibers" 7r~^(p) need not be isomorphic to 
each other and U{J^). This can be seen as follows. A general density matrix is by 
definition a positive operator and can accordingly have null eigenvalues. Consequently, 
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an operator W projected to a given p is not necessarily of full rank. Moreover, if Wi 
and W2 are projected to pi and p2, respectively, they can differ in rank. Hence, in the 
presence of null eigenvalues of p, we can expect a one-to-one correspondence for the 
elements of 7T~^{p) only to a subgroup of U{Jif). Therefore, in order to obtain a well 
defined ?7(^)-bundle with all fibers isomorphic to U{Jif), we need to restrict the base 
space to strictly positive operators (only non-null eigenvalues), 

= {p E M\p> 0}. (16) 

Now we need a subspace S{J^^^) C S{J€'^^) which projects to under vr. We find 
this subspace to be 

S{^^^) = {W e S{jr^^) I KeT{W) = 0}. (17) 

Among others, we just excluded the projective Hilbert space from . But it turns 
out that, once the metric tensor field we are interested in has been derived on Ai^, it 
can be extended to equip the entire quantum state space A4 with a Riemannian metric 

m- 

To obtain the decomposition into horizontal tangent spaces and vertical tangent 
spaces Vw, TwS{,^^^) = Hw ® Vw, we can again introduce a connection to the U{J^) 
bundle S{J^^^) using the scalar product on Hilbert-Schmidt space. Note that 

the tangent spaces Ty/S{^^^) at a point W G S{J^^^) can be identified with subspaces 
of e^"^ due to the Hilbert -space property of Jif^^. For the push forward of a vector 
X e TwS{JF^^) we obtain 

vr,(X) = WX^ + XW^ e Tww^M+. (18) 

Since vr eliminates all the vertical directions, a vector Y is vertical if 

WY^ + YW^ = 0. (19) 

For horizontality of X, we thus require (X, F)hs = to hold for all vertical vectors 
Y E Vw- This leads to the condition 

X^W - WX^ = 0. (20) 

Note that X, being a tangent vector to the point W, can equivalently be written as 
^W{t)\t=o for some curve W{t), W{0) = W. Then, following [211 122], one can show that 



the ansatz dW = GW, with G a hermitian matrix-valued 1-form, solves equation (20). 
For the 1-form dp, defined as the exterior derivative of the density matrix p = WW\ 
this translates into 

dp = pG + Gp. (21) 

If we now define a metric tensor field on the tangent spaces of by taking, again, 
the real part of the scalar product and admitting only horizontal vectors as arguments, 
we obtain the so-called Bures metric 

g{Pi, P2) = m^i, X2)bs} = dp ® G{Pi, P2), (22) 
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with Pi = 7r*(Xf^), P2 = 7r*(X|^). Solving equation (21) for the matrix elements of G 
and using a spectral resolution of the identity operator 1 = J2n \'^n){4'n 
eigenvectors of the density operator p, one obtains 

2 ^ 



9 



in terms of the 



(23) 



Pn+Pm 

for the metric tensor field, which was first found by Hiibner [23]. The Pn are the 
eigenvalues of the density operator p, which can be interpreted as statistical weights. In 
reference [21], by expanding dp in terms of the eigenstates of p, equation (23) is taken 
as a starting point for decomposing g into two parts, 

(24) 



9 



with 



and 



9' 



4 ^ 

n 

-Ay. 



'Pn 



-Apn 



-Apn 



jPn - Pmf 
Pn ~l~ Pm 



iPr, 



(25) 



(26) 



The so-called classical (cl) contribution g^"^ formally coincides with the Fisher-Rao metric 
of classical information geometry [24j. g^'^ ^ in contrast, was dubbed the non-classical 
(nc) contribution. In reference [21] it was also shown that the Bures metric indeed 
reduces to the Fubini-Study metric for pure states. 

As a last step, we need to argue that the Bures metric defined on can be 
extended to M.: An explicit calculation for finite systems reveals that the subspaces 
corresponding to zero-eigenvalues do not contribute to the trace operation which finally 
yields the distance between two density operators. Hence, equations (25) and (26) 
can be continued to M. without modifications, and we have successfully constructed a 
fidelity metric on the quantum state space M.. The two expressions (25) and (26) form 
the starting point for our discussion of phase transitions at finite temperature and their 
relation to the Riemannian structure of quantum state space. 



4. Isotropic LMG Model 

In this section we introduce the Lipkin-Meshkov- Glick (LMG) model in its spin 
formulation and give some of its basic properties, following mainly the presentation 
in [2S]. We then specialize to the isotropic case which is exactly solvable with little 
effort even in the case of finite systems and shortly report on its ground-state structure. 
Finally, the exact thermodynamic solution is recalled. 

The LMG model describes spin-1/2 degrees of freedom residing on the vertices 
of a graph. The spins interact through a ferromagnetic exchange coupling of infinite 
range, i.e., all spin pairs interact with equal strength, 

^LMG = E ^^^""^ - ^ E (27) 
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where 

cT. = ia:,ala:)\ <t\ = [a^ , at) . (28) 

Here, i,j label the graph vertices and cTj denotes the vector of Pauli matrices acting 
on the Hilbert subspace = corresponding to each vertex. The Pauli-vector 
components satisfy [o"f,(Tj] = 2i 5 ije ^^a'^ , where we used Greek indices to label spatial 
vector components. The full Hilbert space is given by the tensor product 

N 

= (0 J^i = [CY^ . (29) 

i=l 



The coupling matrix C in (27) is given by C = diag(l,7,0), with anisotropy parameter 



7. A factor of 1/N is included in (27) to ensure a finite free energy per degree of freedom 
when taking the thermodynamic limit. Moreover, an external magnetic field of strength 
h, pointing in the z-direction, tries to align the spins along this direction. 

The model dynamics can be formulated entirely in terms of the total spin 

1 ^ 

S=-J2^^- (30) 

1=1 

Its components obey the usual angular momentum commutation relations [S*^, S^] = 
ie^^^S^, yielding [S"^, S^] = for all /x. Introducing spin-raising and -lowering operators 
S± = {Sx ± iSy) /2, the Hamiltonian can be rewritten as 

^LMG = (s^ - Si - I) - 2hSz -^-^{Sl + si) . (31) 

Since [S^, TYlmg] = 0, is a conserved quantity under the dynamics induced by 7-^lmg- 
Thus, the Hilbert space can be decomposed as 

^={€'f'' = @dsVs. (32) 

s 

where ds denote the multiplicities of irreducible and unitary S'f/(2)-representations Vs 
of dimension dim 1)5 = 25* + 1. For convenience, we choose even in the following, 
obtaining S G {0, . . . , N/2}. Moreover, the Hamiltonian is invariant under time reversal 
{h I— >■ —h,cr I— i> — cr). Therefore, all eigenvalues are at least twice degenerate (Kramers 
degeneracy), En{h) = En'{—h), where n and n' denote distinct sets of quantum numbers. 
Due to this symmetry, we can restrict the discussion of the spectral properties of T^lmg 
to the case h ^ 0. We now specialize to the isotropic model, and comment on the 
anisotropic case in section [7j The Hamiltonian (31) reduces in the isotropic case 7 = 1 
to 

ntuG = -If (s' - si - ^ V '^hSz. (33) 



N \ ^2 



Since [5'2,7iLMG] — 0' ^^"^ ^^^o Sz is an integral of motion. We denote by \SM) the 
simultaneous eigenstates of 5*^ and 5*2, where 

S'^\SM) =S{S ^\)\SM), Sz\SM) = M\SM). (34) 
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For every spin sector Vs, the angular momentum eigenstates \SM), M = —S, . . . , +5, 
are eigenstates of the Hamiltonian, and T^lmg invariant under rotations about the 
z-axis. Its eigenvalues are given by 

= -^J S{S + 1) - ~2hM. (35) 



EsM attains its minimum in the maximum-spin sector (i.e., for quantum number 
5*0 = A^/2) with magnetic quantum number 

rx(.iV/2) fa 

I N/2 for /i ^ 1, 

where 

( \x\ for X = H +5, 5 G [0,1/2), 

[\x] for x=[x\+S,Se [1/2, 1), 

is the rounding function. The value of the external field h therefore determines which 
of the angular momentum states G is selected as the ground state. At certain 

values of h the ground state switches from one M- value to another (see equation 36), 
and these points of degeneracy are termed level crossings. In the thermodynamic limit 
N ^ oo, the ground-state energy per spin converges towards a continuous function 
of h [25j, being infinitely differentiable almost everywhere. Only at h = hc\T=o = il 
its second derivative with respect to h is discontinuous, signaling the above-mentioned 
phase transitions. 

We now recall the exact thermodynamic solution of the LMG model, which is a 
special case of a result by Pearce and Thompson [2S] obtained for a large class of mean- 
field type spin models in an external field. For the isotropic LMG model, the free energy 
per spin in the thermodynamic limit A^ — oo is given by 

fi(3,h) ^ - lun^^lntTexpi-mtuG) = - r Mn (2 cosh [(3 ^ fi% + h^)) , (38) 

where f3 denotes inverse temperature. The relative magnetization in z-direction, 
fiz = —df/dh, is completely determined by the value of the external field h and the 
scalar order parameter fi^y = fJ^xyiP, h). The latter obeys the self-consistency equation 

+ h'= (tanh + h^) ) (39) 

and has the interpretation of a relative in-plane magnetization with respect to the 
maximum total spin A^/2, 

^^,=jL-]^((i:<)^(i:-.f) («) 

where (■) denotes a thermal-equilibrium average. 

The self-consistency equation (39) determines the phase diagram completely, see 
figure [1] For fields with \h\ < 1 and temperatures T below the critical temperature 
, Hxy takes non-zero values and vanishes continuously when approaching the phase 



Finite- Temperature Fidelity- Metric Approach to the Lipkin-Meshkov- Click Model 10 




Figure 1. Phase diagram of the isotropic LMG model plotted in the (T, /i)-plane. The 
blue shaded area marks the ordered phase with finite in-plane magnetization 7^ 0, 
separated from the paramagnetic phase fx^y = (shown in black) by a line of phase 
transitions making up the phase boundary (dashed white). 



boundary by either an increase in temperature or magnetic field. The phase boundary 
reached as fi^y = consists of all points {f3c, he) that obey he = tanh(/?c^c) or 

Pc{hc) = /i^"'^arctanh(/ic), (41) 

where /3c = 1/7"^ and he denote critical values of inverse temperature and magnetic 
field. In summary, the LMG model shows a ferromagnetically-ordered phase separated 
from a paramagnetic phase by a line of continuous phase transitions. All these exact 
thermodynamic results for the infinite system coincide with results obtained from a 
mean- field treatment as reported in [11] and [27] . 



5. Metric Tensor Field for Thermal States 



In this section we compute the metric tensor field on the submanifold of thermal states. 
A thermal equilibrium state (or Gibbs state) of the isotropic LMG model is given by 

where 

N/2 +S 

Zj^iP, h) = J2ds ^M-H^'-s^-fh^'^s^) \SM) (43) 

5=0 M=-S 

is the canonical partition function. Here, the density operator inherits a dependence on 
(3 and h from the Hamiltonian and the partition function. 

For finite systems, equation (42) defines a parameterization of the submanifold Q 
of thermal states. Equivalently, we can take this as a trivial chart p{f3,h) i— {(3,h), 
defining local coordinates on Q. Vector fields (and, analogously, 1-forms or higher rank 
tensor fields) can then be expressed with respect to the coordinate basis {9/3, dh}- 
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As a first step, we use the decomposition = ^gdsVs of the Hilbert space 
of the LMG model to cast the spectral representation of g [equations (24)-(26)] in a 
different form: We solve equation (21), separately within every spin sector Vs, for the 
matrix elements of G with respect to the (25' + l)-dimensional basis {\SM)}. Plugging 
the result into equation (|22|) and taking the trace, we arrive at the expressions 



7 , ^ dpsM ® , ^ dpsM (44) 

4 ^ ^ y/p^ y/p^ 



and 



g- = lYdsY (Pg^ PSM') ^sM\d\SM') ® {SM\d\SM'y = 0. (45) 

1^' PSM+PSM' 



Here, by 



PSM = ^ exp{-f3EsM) (46) 
An 



we denote the statistical weights with energies Esm as given in (35). Since the 
eigenstates of T^lmg do not carry any explicit /i-dependence (nor, of course, any /?- 
dependence), the non-classical contribution g"^"^ vanishes by virtue of d \ SM) = 0, and 
we are left with the classical Fisher- Rao contribution (44). This is maybe not too 
surprising: Since the operators S"^, in the Hamiltonian of the isotropic LMG model are 
commuting, a (classical) probability distribution can be assigned, and the corresponding 
information geometrical metric is known to be the one of Fisher-Rao. A more detailed 
and more general discussion of Hamiltonians consisting of commuting summands and 
the implications on the fidelity can be found in reference [28]. 

It is straightforward to compute the 1-forms dpsM in equation (44) with respect to 
the dual coordinate basis dj3, dh, yielding 

dpsM = dfiPsAidp + dhPsAidh, (47) 

where the partial derivatives can be rewritten as 

dpPsM = PSM {{ntuG > - Esm) , (48 a) 

dhPsM = -2(3psM ((5,) - M) . (486) 

Expanding the metric tensor field with respect to the rank two tensor basis d/5 ® d/3, 
dp ® dh, etc., we obtain 

*.-iE^»E^^. («") 

(496) 



3hh ^ ^ 

S M 

1 , ■ 

9hl3 



1 \^ , \^ dpPsMdhPSM f.f. N 

zz^^^z^ — ' (49^) 
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and, furthermore, gi3h = Qhp- Inserting (48a) and (486) into the above equations, we 



find the metric components to be given by equihbrium fiuctuations and correlations 
1 
4 

i2 



9pp 

Qhh 



/-T/iso \2 



-T/lSO 

'T-LMG 

Analogous results have been derived in [2l] for the case of a quantum Ising spin chain. 
Since T^lmg contains only commuting operators, the connected correlation functions 



/ISO 

'-LMG 



(50 a) 
(506) 
(50c) 



on the right hand sides of (50a)-(50c) can, up to some pre-factors, be expressed as 
derivatives of the free energy ( [38) ), yielding 

N 

4^^ 

T 

Nf3, 



9PP 



9hh — — 



9hf3 



-dl iPf) , 



(51a) 
(516) 
(51c) 



Note that these expressions suggest a close relationship to Ruppeiner geometry 
[10\ . where a Riemannian metric is defined in terms of derivatives of a suitable 
thermodynamic potential, e.g. the free energy. In the case of a vanishing non-classical 
part g"^^, the Bures metric for thermal states and the Ruppeiner metric indeed differ 
just by a coordinate transformation. 



The exact result for the free energy per spin (38) and the self-consistency condition 



(|40|) can now be used to compute the metric per spin in the thermodynamic limit, 
lim^v- 



9 



(oo) 



N 



g. One would expect from equations (51a)-(51c) that, in this limit 



the metric inherits a nonanalytic behaviour from the nonanalyticity of the free energy 
/, but we will see in the following that the situation is even a bit more involved. 



5. 1 . Ordered phase 

For the ordered phase with 7^ we can collect the metric components in the form 
of a diagonal matrix. 




1 
4 



[^lly+h■^)(l+^lly+h'^ 

l-/3[l-(/.2,+h2)] 





(52) 



Equation (52) provides the components of a well-defined Riemannian metric, and we 



can now interpret the metric components as indicators of how well thermal states with 
close-by values of (3 and h can be distinguished. The graphs of g^^'' and g^^^ in the 
(T, /i)-plane are shown in figure [2] 

Instead of considering the metric components separately, one can combine them to 
compute the Ricci scalar, a quantity characterizing the curvature of a manifold. We 
have computed the Ricci scalar by making use of the Maurer-Cartan p9j equations. 
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Figure 2. Components g^^ and g^^f^ of the metric in the thermodynamic Hmit as 
functions of temperature T and external magnetic field h. When approaching the phase 
boundary, g'j^^ increases, suggesting enhanced distinguishability under variations of 

temperature. From the divergence of gj°^^ at T = one can see that the Bures metric is 
not well-defined for the ground state in the thermodynamic limit (nor is it well-defined 
for finite systems at T = 0) . 



Here, we skip the details of this calculation and present only the final result in the 
thermodynamic limit, 



PQ\ Q P 



(53) 



where P = \/j3/2 and Q = ^J Hxydpl-^xyl^. The graph of 7^*^°°) is shown in figure 3j 
We observe that 7^*^°°^ is negative in the entire ordered phase. In reference [23], it 
was conjectured that a negative Ricci scalar should correspond to the "classical realm" 
of a given system. Since the Hamiltonian of the isotropic LMG model contains only 
mutually commuting terms, it may be regarded as classical, and the negative Ricci 
scalar we observe is in agreement with the conjecture. 

5.2. Paramagnetic phase 

Surprisingly at first sight, the metric not only becomes singular at the phase boundary, 
but changes its structure entirely from one phase to the other. Writing the rescaled 
metric tensor in matrix form, 

^i(=°*H"(S (^^> 

we find that in the disordered phase with = 0, the matrix of the metric components 
has vanishing determinant. This in turn implies that in the disordered phase the rank 
two tensor becomes degenerate and is not a proper Riemannian metric anymore. Since 
we started out with a Riemannian metric in the finite-system case, the limit N ^ oo 
must have destroyed this property. 

One can understand the physical origin of this effect by considering expression 




(38) for the free energy of the LMG model: For vanishing in-plane magnetization 
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Figure 3. Ricci scalar 7?,^°°' of the isotropic LMG model in the thermodynamic 
limit as functions of temperature T and external magnetic field h. As for the metric 
components in figure [2] 71'^°°^ is well-defined only for the ferromagnetically ordered 
phase, and the breakdown of its existence can be interpreted as a signal of the phase 
transition. Note that, where defined, TZ^°°^ is negative, suggesting a classical-type 
behaviour of the system. 

fJ'xy = 0, the free energy is identical to that of a spin system coupled to an external 
field h, but without any spin-spin interaction whatsoever. Such a system is governed 
by the Zeeman-Hamiltonian 7i = —hSz, and the corresponding thermal state is given 
by p = exp{l3hSz) with partition function Z = tr exp(/3/i5'^). For this system, 
all {(3, h) with f3h = const, parametrize the same density operator. The very same 
situation occurs also for the paramagnetic phase with n^y = of the LMG model in the 
thermodynamic limit, leading to the mentioned degeneracy of the matrix of the metric 
components. As a consequence, thermal states of this phase should be parametrized by 
only a single parameter, namely the reduced field h = f3h. The corresponding metric 
on such a one-dimensional manifold can be obtained by a computation similar to the 
two-dimensional case reported above, yielding 

= (4 cosh h) All ® All. (55) 



6. Pubini-Study Limit 

We have mentioned at the end of section |3] that the Bures metric reduces to the 
Fubini-Study metric when considering pure states. This should in principle allow us to 
define a Fubini-Study metric on the ground-state manifold parametrized by a ground- 
state mapping h i— > \iI)gs)- However, we have observed in section |4] that the ground 
state of the isotropic LMG model is the angular momentum eigenstate \N/2, MqIH)) 



where, according to equation (36), Mo(/i) is selected by the rounding function X. As 
a consequence, no different iable parametrization of the ground state exists. Since this 
property naturally carries over to the respective chart mapping, the Fubini-Study metric 
on the ground-state manifold of the isotropic LMG model in the finite system is not 
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well defined. Similar issues, related to ground state level crossings and their effect on 
the fidelity, are also discussed in reference [28j. 

In order to investigate the ground state behaviour in the infinite system, we can 



alternatively study the limit T — >^ of the Riemannian metric characterized by (52). 
The component gj^'' is found to vanish in this limit, in agreement with the fact that 



according to equation ( |50a[ ), it is proportional to the specific heat. The component gl^\ 
however, diverges for T — > asymptotically as T~^. Hence the (rescaled) ground-state 
metric is not well defined in this limit, although the ground state energy becomes a 
continuous function of the external field h in the thermodynamic limit. 

7. Remarks on the Anisotropic Case 

The problems we encountered in the previous section when trying to compute a metric 
on the ground-state manifold of the isotropic LMG model can be traced back to the 
fact that the Hamiltonian consists only of mutually commuting terms, which in turn 
allows for level crossings. For the anisotropic case, in contrast, we would expect a well- 
defined ground state metric. Exact results for the spectrum and eigenvalues, exist also 
for the anisotropic LMG model [301 El]- principle, these results would allow one 
to compute the metric on the manifold of thermal states, but unfortunately, they are 
expressed as rather complicated multiple sums, with coefficients given as solutions of 
differential equations, which makes the calculation quite difficult in practice. A more 
accessible result is given in reference [55] for the energy spectrum of this model in 
the thermodynamic limit but, as far as we can see, this might not be sufficient for a 
computation of the fidelity metric. 

Alternatively, one might try to compute the metric in mean-field approximation. 
Knowing that, in the thermodynamic limit, the mean- field solution of the (isotropic or 
anisotropic) LMG model coincides with the exact solution, one might hope to obtain the 
exact metric from a mean-field calculation as well. Unfortunately this is not the case, 
nor are we aware of any other approximation that retains enough of the original quantum 
state space structure in order to deliver an accurate description of the underlying 
geometry. 



In contrast to the isotropic LMG model where only the classical part (44) was 



found to contribute to the metric, we expect the anisotropic case to have a non- 



zero non-classical contribution (45). Furthermore, it will not anymore be possible to 
completely express the metric in terms of derivatives of the free energy, and only then 
the characterization of phase transitions by means of the fidelity metric would really go 
beyond a thermodynamic description. 

8. Conclusions 



In this paper we have followed the idea that a suitable metric on quantum state space can 
be used to identify and characterize both classical and quantum phase transitions. We 
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have reviewed how such a fidehty metric is constructed, either on the space of pure states 
V{J^) or on the space of state operators Ai. From the Bures metric, i.e., the fidehty 
metric on Ai, the metric tensor field on the submanifold of thermal states has been 
derived. As an application of these concepts, we studied the LMG model of spin-1/2 
degrees of freedom sitting on the vertices of a fully connected graph. The choice of this 
model was mainly motivated by the fact that its thermodynamics is particularly simple 
to solve, and exact results are available for the free energy per spin in the thermodynamic 
limit, both for the isotropic and the anisotropic LMG model. 

For the isotropic LMG model, we computed the metric tensor field on the 
submanifold of thermal states, and we found that all metric components can be written 
as derivatives of the free energy. This implies a close relation to Ruppeiner geometry, 
but this should be a peculiarity of models with purely classical contributions to the 
metric. Another peculiar feature special to the isotropic case is that on the ground- 
state manifold the metric is not well defined, neither by direct construction from the 
finite system nor by a detour via thermal states and the subsequent zero-temperature 
limit. This can be seen as a consequence of level crossings which occur in this case, but 
are avoided in the anisotropic model. 

As expected, we find that the phase transition of the isotropic LMG model occurring 
at the transition line (41) in the (T, h)-plane is well captured by the metric components. 
In a way, this signature is even more pronounced than for other models which had been 
studied before: Not only do the metric components show a singularity or discontinuity, 
but we find that, in the thermodynamic limit, the tensor field on the (T, h)-plane 
becomes degenerate for the paramagnetic phase and therefore ceases to be a proper 
Riemannian metric. 

It would be worthwhile to compare these results to the corresponding metric of the 
anisotropic LMG model. Here, we expect the metric to be well-defined on the ground- 
state manifold, and the non-classical part (45) of the metric to give a non- vanishing 
contribution. 
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